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Abstract 

Spatially homogeneous and anisotropic Bianchi type VIq cosmo- 
logical models with cosmological constant are investigated in the pres- 
ence of anisotropic dark energy. We examine the effects of electromag- 
netic field on the dynamics of the universe and anisotropic behavior 
of dark energy. The law of variation of the mean Hubble parameter 
is used to find exact solutions of the Einstein field equations. We 
find that electromagnetic field promotes anisotropic behavior of dark 
energy which becomes isotropic for future evolution. It is concluded 
that the isotropic behavior of the universe model is seen even in the 
presence of electromagnetic field and anisotropic fluid. 
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1 Introduction 



The most remarkable advancement in cosmology is its observational evidence 
which says that our universe is in an accelerating expansion phase. Supernova 
I a data [U [2] gave the first indication of the accelerated expansion of the uni- 
verse. This was confirmed by the observations of anisotropics in the cosmic 
microwave background (CMB) radiation as seen in the data from satellite 
such as WMAP [3] and large scale structure [I]. Today's one of the major 
concerns of cosmology is the dark energy. Recent cosmological observations 
PQ |2j [33 H] suggest that our universe is (approximately) spatially flat and 
its cosmic inflation is due to the matter field (dark energy) having negative 
pressure (violating energy conditions). The composition of universe density 
is the following: 74% dark energy (DE), 22% dark matter and 4% ordinary 
matter [5]. Though there is a compelling evidence that expansion of the 
universe is accelerating, yet the nature of dark energy has been under con- 
sideration since the last decade 0-110] • Several models have been proposed 
for this purpose, e.g., Chaplygin gas, phantoms, quintessence, cosmological 
constant and dark energy in brane worlds. However, none of these models 
can be regarded as being entirely convincing so far. 

The cosmological constant, A is the most obvious theoretical candidate of 
DE which has the equation of state (EoS) u = —1. Astronomical observations 
indicate that the cosmological constant is many orders of magnitude smaller 
than estimated in modern theories of elementary particles [IT] . Stabell and 
Refsdal [12] discussed the evolution of Friedmann-Lemaitre Robertson and 
Walker (FLRW) dust models in the presence of positive cosmological con- 
stant. These results are given in a more generalized form [T3T [T4] using the 
general EoS. Wald [T5] examined the late time behavior of expanding homo- 
geneous cosmological models satisfying the Einstein field equations (EFEs) 
with a positive cosmological constant. He found that all the Bianchi type 
models except IX showed isotropic behavior. These models exponentially 
evolve towards the de Sitter universe with a scale factor (3/A) 1 / 2 . Goliath 
and Ellis JTB] used the dynamical system methods to observe the spatially 
homogeneous cosmological models with a cosmological constant. The in- 
clusion of cosmological constant provides an effective mean of isotropizing 
homogeneous universes [151 US] • 

The presence of magnetic fields in galatic and intergalatic spaces is ev- 
ident from recent observations [T7]. The large scale magnetic fields can be 
detected by observing their effects on the CMB radiation. These fields would 
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enhance anisotropies in the CMB, since the expansion rate will be different 
depending on the directions of the field lines [TSJ [TH] - Matravers and Tsagas 
[20] found that interaction of the cosmological magnetic field with the space- 
time geometry could affect the expansion of the universe. If the curvature 
is strong, then even the weak magnetic field will effect the evolution of the 
universe. The magneto-curvature coupling tends to accelerate the positively 
curved regions while it decelerates the negatively curved regions [201 121] • 

Jacobs [22] studied the spatially homogeneous and anisotropic Bianchi 
type I cosmological model with expansion and shear but without rotation. 
He discussed anisotropy in the temperature of CMB and expansion both with 
and without magnetic field. It was concluded that the primordial magnetic 
field produced large expansion anisotropies during the radiation-dominated 
phase but it had negligible effect during the dust-dominated phase. Dunn 
and Tupper [23] discussed properties of Bianchi type VIq models with perfect 
fluid and magnetic field. Roy et al. [21] explored the effects of cosmological 
constant in Bianchi type / and VIq models with perfect fluid and homoge- 
neous magnetic field in the axial direction. They found that model expanded 
for negative values of the cosmological constant while it contracted for pos- 
itive A. In a recent paper [25], Sharif and Shamir explored the vacuum 
solution of Bianchi I and V models in f(R) gravity. 

Rodrigues [26] proposed Bianchi / with a non-dynamical DE component 
which yields anisotropic vacuum pressure in two ways: (i) by considering the 
anisotropic vacuum consistent with energy-momentum conservation; (ii) by 
implementing a Poisson structure deformations between canonical momenta 
such that re-scaling of scale factors is not violated. Koivisto and Mota [2Z] 
have investigated a cosmological model containing the DE fluid with non- 
dynamical anisotropic EoS and interacts with perfect fluid. They suggested 
that if the DE EoS is anisotropic, the expansion rate of the universe becomes 
direction dependent at late times and cosmological models with anisotropic 
EoS can explain some of the observed anomalies in CMB. 

Recently, Akarsu and Kilinc [2H] investigated anisotropic Bianchi type / 
models in the presence of perfect fluid and minimally interacting DE with 
anisotropic EoS parameter. They found that anisotropy of the DE did not 
always promote anisotropy of the expansion. The anisotropic fluid may sup- 
port isotropization of the expansion for relatively earlier times in the universe. 
The same authors [29J have worked on the Bianchi type III model in the 
presence of single imperfect fluid with dynamical anisotropic EoS parameter 
and dynamical energy density. They observed that anisotropy of the expan- 
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sion vanished and hence the universe approached isotropy for late times of 
the universe in accelerating models. 

It would be worthwhile to see what happens if we consider the electro- 
magnetic field with anisotropic DE. The main purpose of this work is to 
look at the effects of electromagnetic field on the dynamics of the universe 
in the presence of anisotropic DE for Bianchi type VIq. The layout of the 
paper is as follows. In section 2 we describe spatially homogeneous and 
anisotropic Bianchi type VIq spacetime and formulate the EFEs in the pres- 
ence of anisotropic fluid and magnetic field. Section 3 presents a special 
law of variation for the mean Hubble parameter which yields constant de- 
celeration parameter. This law generates two types of solutions: power law 
and exponential expansion. In section 4, a hypothetical form of fluid is ob- 
tained by making an assumption on anisotropy of the fluid. We obtain exact 
solutions of the EFEs and discuss physical behavior of anisotropic DE and 
universe model. Finally, section 5 concludes the results. 



2 Model and the Field Equations 

The spatially homogeneous and anisotropic Bianchi type VIq model is de- 
scribed by the line element 

ds 2 = dt 2 - A 2 (t)dx 2 - e 2mx B 2 (t)dy 2 - e~ 2mx C 2 (t)dz 2 , (1) 

where scale factors A, B and C are functions of cosmic time t only, m ^ 
is a constant. The energy- momentum tensor for the electromagnetic field is 
given as [30J 

T^=fc[\h\ 2 (uX-l/25;)-h,h% (2) 
where u v is the four-velocity vector satisfying 

g^u»u v = 1. (3) 

Jt is the magnetic permeability and is the four-magnetic flux given by 

K = ^fe^F^u", 0, v,a,P = 0, 1, 2, 3) (4) 

where e^ap is the Levi-Civita tensor, F Q/3 is the electromagnetic field tensor 
and \h\ 2 = h u h u . We assume that magnetic field is due to an electric current 
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produced along x-axis and thus it is in yz-plane. In co-moving coordinates 
u v = (1,0,0,0) and hence Eq.(4) gives hi ^ 0, h = h 2 = h 3 = 0. Using 
these values in Eq. (SJ), it follows that F\ 2 = Fi 3 = 0, F 23 ^ 0. 

The electric and magnetic field in terms of field tensor are defined as 



Efj, = F^u", = -e lluol F va . (5) 

According to Ohm's law, we have 

h^ v J u = oF ilv u h ', (6) 

where h^ v = g^ u + u^u u is the projection tensor orthogonal to u^, a is the 
conductivity and J M is the four current density. In the magnetohydrodynamic 
limit, conductivity takes infinitely large value while current remains finite so 
that Efj, — y [31]. Consequently, Eq.(|5]) leads to F i = F 2 = F 03 = 0. Thus 
the only non- vanishing component of electromagnetic field tensor F^ is F 23 . 
The Maxwell's equations 

Ffj,i/-a F va -^ -\- F a ^- V 0, (7) 

are satisfied by 

F23 = K = constant. (8) 

It follows from Eq.(j4]) that 

AK ,„ K 2 



Using this equation in Eq.(j2J), we obtain 

K 2 
2J[B 2 C 

Thus we have 



rpO(ern) ^ rpl(em) rp2(em) rp3(em) 



K 2 K 2 K 2 K 2 

T u(em) = diaQ \_^ * * * 1 QO) 

" 9[ 2jlB 2 C 2 ' 2jlB 2 C 2 ' 2JlB 2 C 2 ' 2jlB 2 C 21 ' { } 

The energy-momentum tensor for anisotropic DE fluid is taken in the 
following form [28| |29] 

T v = diag[p, -p x , -p y , -p z ] (11) 
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This model of the DE is characterized by the EoS, p = cop, where co is not 
necessarily constant |33J. From Eq.f llip . we have 



T; = diag[l, -(co + 5), -co, -(co + 7 )]p, (12) 

where p is the energy density of the fluid; p x , p y and p z are pressures and 
co x , co y , co z are directional EoS parameters on x, y and z axes respectively 
The deviation from isotropy is obtained by setting 

co x = co + 5, co y = co, co z = co + 7, 

where co is the deviation free EoS parameter and 5 and 7 are the deviations 
from co on x and z axes respectively. The EFEs with cosmological constant 
are given by 

Gfj, v = Rfj,u - \R9^ ~ = 8vr(T MV + T^ m) ), (13) 

where R^ v is the Ricci tensor, R is the Ricci scalar, T^ u is the energy- 
momentum tensor for anisotropic fluid and T^J^ is the energy-momentum 
tensor for the electromagnetic field. 

For Bianchi type VIq spacetime, the EFEs become 

AB AC BC m 2 AttK 2 . , , 

AB + AC + BC-A^ = MP+ WC 2+K (M) 

B C BC m 2 . AttK 2 a , . 

B + C + BC + ir 2 =^ + 5)p+ Wc 2+ ^ (15) 

A C AC m 2 o AttK 2 . 

A + C + AC-^ = -^ UJp -Wc 2+K (16) 

A B AB m 2 , . AttK 2 4 , . 

A + B + AB-^ = - 8 ^ + -> :) "-W^ + A ' <17) 

m(§-§) = 0, (18) 
where dot denotes derivative with respect to time t. Equation (fl8|) yields 

B = Cl C, (19) 



where c% is a constant of integration. Subtracting Eq.l fP7|) from ( Tl6|) and 
using ( fl9l) . we obtain 7 = which shows that directional EoS parameters 
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Uy, u z and the pressures p y , p z become equal. Using Eq. (IT5|) and 7 = 0, the 
EFEs (IT4"|) - (IT7)) reduce to the following set of equations 

AB B 2 m 2 AtxK 2 , , . 

2 AB + W-X = 8 " >+ pfc?F + A - (20) 

B B 2 m 2 , rN 47ri^ 2 , , s 

2 B + W + X = - 8 ^ + i »" + =jBl + A > < 21) 

A B AB m 2 AttK 2 t . , 

A + B + AB-lV = ~^"-pW + A - <22) 

3 Some Physical and Geometrical Parame- 
ters 

Here we discuss some physical and geometrical quantities for the Bianchi type 
VIq model which are important in cosmological observations. The average 
scale factor is given by 

a = {k x AB 2 )^ (23) 



while the volume is defined as 

V = a 3 = k x AB 2 . (24) 

The mean Hubble parameter H and the directional Hubble parameters Hi(i = 
1, 2, 3) in x, y and z directions are 

ff=I(lnVj=lna = i(^ + 2|), H x = j, Hy = H z = ^. (25) 

The physical parameters such as scalar expansion O, shear scalar a 2 and 
anisotropy of the expansion A are given as follows 

6 = u% (26) 
a 2 = l -a ah o a \ (27) 

3 

I 

A 



^E(^V^ )2 - (28) 

The anisotropy of expansion shows isotropic behavior for A = 0. 
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It is mentioned here that any universe model becomes isotropic for the 
diagonal energy-momentum tensor when t — > +00, A — > 0, V — > +00 and 
T 00 > (p > 0) [291 E2]- The law of variation of mean Hubble parameter is 
given as 

H = la- n = IthAB 2 )-^ 3 , (29) 

where I > and n ^ 0. This law was initially proposed by Berman [M] 
for spatially homogeneous and isotropic RW spacetime which yields constant 
value of the deceleration parameter. In recent papers [28l 135] . a similar law 
is proposed for the homogeneous and anisotropic Bianchi models to generate 
exact solutions. 

The volumetric deceleration parameter q is the measure of rate at which 
expansion of the universe slows down due to self-gravitation. It is defined as 

ad , , 

" - - 4T < 30 > 

Using Eqs.(E>D and (EHD, we get 

a = la' n+1 , a = -l 2 {n - l)aT 2n+l (31) 

which gives constant values of the deceleration parameter as follows 

q = n — 1 for n^O, 
q — — 1 for n = 0. 

Recent observations show that expansion rate of the universe is accelerating 
which may be due to the presence of DE and q < 0. Thus the sign of 
q indicates whether the cosmological model inflates or not. For q > (i.e. 
n > 1), the model represents decelerating universe where as the negative sign 
— l^g<0 for 0^72<1 indicates inflation and q = for n — 1 corresponds to 
expansion with constant velocity. Equations fT25|) and fT29|) yield two different 
volumetric expansion laws 

V = k 2 e m , n = 0, (32) 

V = {nlt + c) 3/n , n^0 (33) 

which are used to find exact solutions of the EFEs. In fact these represent 
two different models of the universe. 
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4 Solution of the Field Equations 

We can find the most general form of anisotropy parameter for the expansion 
of Bianchi type VIq in the presence of anisotropic fluid and electromagnetic 
field using Eq.( j25l) . The anisotropy parameter of expansion can be written 

as 

^=^p{H x -H y )\ (34) 

where H x — H y is the difference between the expansion rates on x and y axes 
which can be found using the field equations. 

Subtracting Eq.( f22|) from (|2T|) and after some manipulation, it follows 
that 

A B d I f , c 2m 2 8ttK 2 . , 

A-B-V + Vj iMp+ ^-W^ )Vdt ' <35> 
where d is another constant of integration. Now using this equation in 
Eq. ( I34p . the anisotropy parameter takes the form 



The anisotropy parameter in the presence of isotropic fluid can be obtained 
by choosing 5 = which yields 



2 r , f ,2m 2 8nK 2 , Tr , 2 
A =9IP ld + J^-^ )Vdt]V (37) 

We take the value of 5 so that the integrand in the above equation vanishes 

m 2 K 2 
6 = ~A^A 2+ W^- (38) 

The corresponding energy-momentum tensor for anisotropic DE fluid turns 
out to be 

m 2 K 2 

T; = diagl l, -(„ - — 2 + -u]p. (39) 

The anisotropy parameter of the expansion reduces to 

2 d 2 
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We see that A obtained for the Bianchi type VIq in the presence of anisotropic 
fluid with electromagnetic field is equivalent to that found for the Bianchi 
type III in the presence of anisotropic fluid [29] . 

The difference between the directional Hubble parameters becomes 

H x -H y = ^. (41) 



The most general form of the energy density is found by using Eqs. (1201) and 

(BHD as 

1 r 9/ ^ m 2 4wK 2 . , . 

" = ^ 3H ^-^-^- A -^ ] - (42) 

This shows that anisotropy of expansion, cosmological constant and electro- 
magnetic field reduce the energy density p of anisotropic DE. 

4.1 Model for n = (q = -1) 

The spatial volume of the universe for this model is given by 

V = k 2 e m . (43) 

Using this value of V in Eq.flU) and then solving the EFEs fl20>(p2"j) . the 
scale factors become 

A = k 3 e lt -*—2 e , (44) 

£? = (^-^— ) 2 e 9ife 2 , (45) 
fcifc 3 

~ A&2N.I it+I d e -3K 

(7 = (— — Pe 9ife 2 , (46) 
h 

where k\, k 2 and k 3 are constants of integration. 

The directional and the mean Hubble parameters will become 

H x = l + ~^e- 3 \ H y = H z = l-\U-*\ H = l (47) 



while the anisotropy parameter of the expansion takes the form 

2 d 2 
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A = 5/% e "" (48) 



The expansion and shear scalar are found as 



9 = «• =^ + 2| = 3/ = 3^, (49) 



l r .A.o .,B>„ l„o Id 2 

~2 K 



a 2 = _[(ll) 2 + 2(-) 2 ] - -6 2 = -^e~ m . (50) 
2 lK A J K B J J 6 3k 2 



The energy density of the DE is evaluated by using Eq. (l2"Ui) with the scale 
factors as 



1 r o , Id 2 m m 2 _2Zt+^^- e - 3 " AnK 2 ,k 2 . 
8tt 3 k 2 k 2 -p y k\ 



x e" 41I "^ e ]. (51) 

The deviation free part of anisotropic EoS parameter u can be obtained by 
using Eqs.dHD-(|l6D and (Ell) in Eq. (121)1 



rn; 2,2,2 . r2,2 -6Zt qa;2,2 Q 2,2 -2ZH-f yp-e" 3 " ^VF/Cgf^ 2 

u = {9k 2 k 3 l + a k 3 e — SAk 2 k 3 — 3m /c 2 e 9 ^2 H = 

- 4i( -^ e_3 "l /rj2 7„2^-6tt , qa,Jl2 n/„2,„2,2 , o_27„2 



x e — 9i^ e }/{d 2 ^ e -^« + 3AJfc|fc| - 9k 2 k 2 3 l 2 + 3m 2 k, 



2 



- 2 Zt+4 d e -3M , Yl-Kk%K l _ 4H _4 



d --3(t, 



x e —^9^ + j*" e"" r ~^ e }. (52) 

A 4 



Using the scale factors and the energy density in Eq. (13811. the deviation in 
EoS parameter along x-axis 5 is given as 

, r R 2,2 -2lt+±jfe- 3lt . 1^K 2 k\ _ 4Zt * <L e -3i* 

= {om /c 2 e 9ifc 2 -| — — -e 9 ife 2 }/{« k 3 e + 3Ak| 

7 2 07272,2 ,0 2,2 -2«t+| 7 |- e - 3 » . 127T/c|i^ 2 _4j t _4 <L e -3It 

x k^ — Qk^k^l+Smk^e 9 "=2 -| -f — e 9ife 2 }. (53) 

4.2 Some Physical Aspects of the Model 

We find that the directional Hubble parameters are dynamical where as the 
mean Hubble parameter is constant. Also, the directional Hubble parameters 
become constant at t — and when t — > 00. These deviate from the mean 
Hubble parameter by some constant factor at t — but coincide when t — > 00. 
Since the constant is positive (negative), it increases (decreases) expansion on 
the x-axis and it decreases (increases) expansion on y and z axes. The volume 
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Figure 1: Plot of (a) p(t) and (b) oj(t) verses cosmic time t for m = 1 and 
varying values of I as follows: solid, I = 4; dotted, / = 4.5; dashed, / = 5. 

V of the universe is finite at t — 0, expands exponentially with the increase in 
time t and takes infinitely large value as t — > oo. Thus the universe evolves 
with constant volume and expands exponentially The expansion scalar is 
constant for 0<t<oo and hence the model represents uniform expansion. 

It is mentioned here that the scale factors A(t), B(t) and C(t) are finite 
at t = which implies that the model has no initial singularity whereas 
these diverge for later times of the universe. The anisotropy parameter of 
the expansion and shear scalar are found to be finite for earlier times of the 
universe where as these decrease with time and become zero as t — > oo. This 
shows that anisotropy of the expansion is not supported by the anisotropic 
DE and electromagnetic field. The quantities p, oj and S are dynamical and 
are finite at t — 0. The deviation free EoS parameter of the DE may begin 
in the phantom (oj < —1) or quintessence region (oj > —1) but a;—)- — 1 for 
later times of the universe. One can observe that p increases when u is in 
phantom region and attains a constant value as u— > — 1. This is shown in 
Figure 1. 

We see from Eq. (!53|) that electromagnetic field favors the deviation from 
oj on x-axis, i.e., it contributes to anisotropic behavior of the fluid. The 
deviation parameter increases from negative values towards zero with the 
increase in time and tends to zero as t — » oo shown in Figure 2. Thus the 
anisotropic fluid becomes isotropic for the later times of the universe in the 
case of exponential volumetric expansion. It follows from Eq. fl5Tj) that when 
t — )■ oo, p — > 3l 2 — A = 3H 2 — A which implies that for t — > oo, p > only 
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tf(t) - 4 



0.2 0.4 0.6 0.8 1.0 



Figure 2: Plot of 5 verses cosmic time t for m — 1 and varying values of / as 
follows: solid, / = 4; dotted, / = 4.5; dashed, I = 5. 



if H > a/A/3 and hence A— s-0 and V — > oo. Thus the model approaches to 
isotropy for its future evolution. Here q = —1, dH/dt = which gives the 
largest value of the Hubble parameter and accelerating rate of expansion. 

4.3 Model for n ^ (q = n-l) 

The initial time of the universe is found by using Eq.(l33l) 

£* = -c/nl for n^0. (54) 

We re-define the cosmic time as 

t' = nlt + c (55) 

such that the initial time turns out to be zero, i.e., t' = 0. For this value of 
cosmic time, we can re-define the Bianchi model in the form 

ds 2 = (nl)- 2 dt' 2 - A 2 (t')dx 2 - e 2mx B 2 (t')dy 2 - e - 2mx C 2 (t')dz 2 . (56) 

The corresponding volume will become 



V = t' 



(57) 
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Using this value of V in Eq.QHJ and then solving the EFEs (120)1 - (122]) . the 
scale factors turn out to be 

A{t') = M'-ei^*' 1 " 1 , (58) 

. . . , fe .1 ,1 1 nd j'l-jl 

C(t) = (— )2f«e 3n-3 C . (60) 

^3 

The directional and mean Hubble parameters take the form 

H=(nt')-\ Hx = — + —t'-i, H v = H z = — --t'-^. (61) 
V ; ' nt' 3 y nt' 3 V ; 

The corresponding anisotropy parameter of the expansion turns out to be 

A = \n 2 dH' 2 "t (62) 
9 

while the expansion and shear scalar are 

e = A = 3 # ; <7 2 =ldY-i (63) 
nt 3 

The energy density can be found from Eq. (f20|) by using the scale factors 
Eq.§HD-dnnD as 

... 1 r „. , s 9 , 1 , 9 , 6 m 2 ,2 4 nd ,/1-f 47Ti^ 2 

pt' = 5- 3 <T 2 -A--d¥-" -vji'^e ra* _ __ 

87T 3 K3 /i 

,kl, , 4 4^<l + /i-# 



x Qy-S e tt* x -*]. (64) 



Using this equation and the scale factors in Eq. (12 1 [) . we obtain the deviation 
free EoS parameter u 



w(t') = {(9K)" 2 + dH'-- - 6(nt 2 )- 1 - 3A)A; 2 A; 2 + 3t'~^e^ 



^t^-n 

n — 6 



A A nd ,11 



x ( — J m 2 fc 2 t'ne-3— * ")}/{(dY" + 3A-9(n£')~ 2 ) 

/i 

,9,9 / 4 4 nd j/l-j- ,AnkjK 2 , ,2 8 nd j- , 

x fc 2 A;f + 3t'-« e ra* " ( — J — + m 2 fc 2 t'"e"3— * ")}. (65) 

p 
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Figure 3: Evolution of (a) A(f) and (b) B(f) (dashed line), C(f) (dotted 
line) for n > 3, d > as i! — > oo. 

Finally, the deviation parameter 5 can be obtained using the value of p(t') 
along with Eqs.d58jl-(j60|l in ([38]) 

„, ,. r , 4 4 _nd^,rl-%- AT[kiK 2 9 , 9 ,2 8j!i t 'l-|,, /r , ,9 , 6 

(J(t') = {6t'-W e 3^* "( -1 +m 2 ^t'n e -3^* n )}/{(rf 2 t'-« 

/i 

,, 9, 9 9 ,4 4 nd ATtkiK 2 9,9 ,2 

+ 3A-9(n*0 )^3 + 3* «era* "( — J — +m 2 ^i'" 

x e _ ra* ")}. (66) 



4.4 Some Physical Aspects of the Model 

The universe model accelerates for < n < 1, decelerates for n > 1 and 
expands with constant velocity for n = 1. The mean Hubble parameter, 
shear scalar and directional Hubble parameters are infinite at the initial 
epoch and tend to zero for later times of the universe. If n > 3, d > then 
A(t') takes infinitely large value while both B(t') and C(t') vanish as if — > oo. 
This indicates that spacetime exhibits "pancake" type singularity which is 
shown in Figure 3. If n > 3, d < 0, then A(t') decreases to zero where as 
both B{t') and C{t') continue to increase as t' — > oo, leading to a "cigar" 
singularity shown in Figure 4. For n < 3, the scale factors A(t'), B(t') and 
C(t') become infinite as t' — > oo given in Figure 5. 

We note that spatial volume V is zero at t' — and it takes infinitely 
large value as t' — > oo. The expansion scalar is infinite at f — and de- 
creases with the increase in cosmic time. Thus the universe starts evolving 
with zero volume at the initial epoch with an infinite rate of expansion and 
expansion rate slows down for the later times of the universe. The dynamics 
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Figure 4: Evolution of (a) A{t') and (b) B(f) (dashed line), C{t') (dotted 
line) for n > 3, d < as i! — > oo. 



A(t'), B(t'), C(t') 




Figure 5: Evolution of A(t'), solid line; B(t') dashed; C(t') dotted for n < 3 
as t' — > oo. 



16 



A(t') 



1 . o 

. 8 
. 6 
. 4 
0.2 







2.0 



1.5 



A(t') 1 • o 




0.5 



2 3 
t' 



(a) 



0.0 0.2 0.4 0.6 0.8 
t' 

(b) 



Figure 6: Plot of A verses cosmic time t'; solid line for n < 3 and dotted for 
n > 3 as (a) if -> oo (b) if ->■ 



of anisotropy parameter of the expansion A depends on the value of n. A 
tends to zero as if — > oo and diverges as t' — > for n < 3 and vice versa for 
n > 3 while it remains constant for n = 3 shown in Figure 6. 

Now we determine such values of n which satisfy p > and are suitable 
for the evolution of the universe. We consider magnetic field, anisotropy of 
expansion and cosmological constant for the behavior of energy density of 
the DE. For this purpose, we discuss the following two cases, i.e., n < 3 and 
n > 3. When n < 3, we have p < as t' — > and hence this model is not 
suitable for representing the relatively earlier times of the universe. When 
if — y oo, it follows that p ^ with A ^ which is an appropriate model for 
representing the later times of the universe. For n > 3, we obtain p < as 
t' — > oo and it becomes positive as t' — > 0. Thus this model can represent 
the universe only for the earlier times of the universe by assigning suitable 
values to the constants. 

If n < 3 and 0, p^0 (for later times of the universe), u begins in 
quintessence region and then passes into the phantom region. It remains in 
the phantom region with the increase in time and becomes — 1 as t' — > oo. If 
n > 3 with (for earlier times of the universe), u begins in quintessence 
region, passes into the phantom region for small interval of time then it passes 
back into the quintessence region and becomes — 1 as if — > oo. 

Thus we may examine the behavior of anisotropy of the fluid and universe 
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model for n < 3 with and p^O as t' — > oo. One can observe that mag- 
netic field may increase anisotropic behavior of the DE since it contributes to 
S which decreases with the increase in time and it tends to zero as t' — > oo. 
Thus the anisotropy of the DE vanishes in the presence of magnetic field for 
later times of the universe with negative cosmological constant. We note that 
anisotropy parameter of the expansion is not supported by anisotropy of the 
DE and magnetic field for later times of the universe since for n < 3, A— >0 
as t' — > oo. It is observed that V— >oo and p^O as i! — > oo for n < 3 and 
hence the model represents isotropic universe for future evolution. 

5 Summary and Conclusion 

We have obtained two exact solutions of the dynamical equations for the spa- 
tially homogeneous and anisotropic Bianchi type VIq model with magnetic 
field, anisotropic DE and cosmological constant. The dark energy component 
is dynamical which yield anisotropic pressure. Assuming the law of variation 
of the mean Hubble parameter, the cosmological models are given for n = 
and n^O. The physical and geometrical properties of the models are dis- 
cussed. We have found the explicit form of scale factors and have explained 
the nature of singularities. 

The model represents uniform expansion for the exponential expansion 
while in the case of power law expansion, universe expands with an infinite 
rate of expansion which slows down for the later times of the universe. It is 
found that electromagnetic field affects anisotropies in the CMB, in partic- 
ular, it increases anisotropic behavior. Our results show that even the fluid 
is anisotropic which yields anisotropic EoS parameter with the electromag- 
netic field, its anisotropy vanishes for future evolution of the universe in both 
cases. The expansion of the universe becomes isotropic due to the isotropic 
behavior of the fluid when t' — > oo. 

The model with zero deceleration parameter can approach to isotropy as 
t' — >■ oo with the condition that H > a/A/3, A^O. It is shown that ui is in 
the phantom region which tends to constant value —1 for later times of the 
universe. Thus the expanding model is accelerating. Bianchi models usually 
isotropize for the positive cosmological constant but we have shown that 
model for the power law expansion isotropizes for n < 3, A^O for later times 
of the universe. The above analysis shows that inclusion of the cosmological 
constant in homogeneous cosmological model greatly affects the late time 
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behavior as pointed out by Ellis [IB]- It is interesting to mention here that 
though uj begins in quintessence region then passes into the phantom region, 
but it remains in the phantom region and tends to — 1 as if — > oo which 
can result in accelerating the expansion. However, due to the presence of 
negative cosmological constant, it starts contracting. One can observe that 
this behavior may re-collapse the universe and hence the model represents 
decelerating expansion of the universe. 

Finally, we would like to mention here that the model (the exponential 
expansion law) represents expanding universe for both present and future 
evolution which fits with the current observations and hence the standard 
ACDM model. Also, the observational data favors the power law ACDM 
model pQ . We have found that the model represents an accelerating universe 
for the power law expansion. Both the expansion models approach to the 
EoS of cosmological constant for future evolution. The anisotropy of the 
universe and DE vanish for the period of the accelerated expansion. Thus 
the isotropy is observed for the future evolution of the universe. However, 
there is still a possibility of DE component with anisotropic EoS in present 
epoch of the universe. 
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